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Abstract 

A multiband CuO Hubbard model is studied which incorporates long-range 
(LR) repulsive Coulomb interactions. In the atomic limit, it is shown that a 
charge-transfer from copper to oxygen ions occurs as the strength of the LR 
interaction is increased. The regime of phase separation becomes unstable, and 
is replaced by a uniform state with doubly occupied oxygens. As the holes 
become mobile a superfluid condensate is formed, as suggested by a numerical 
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analysis of pairing correlation functions and flux quantization. Although most 
of the calculations are carried out on one dimensional chains, it is argued that 
the results are also applicable to two dimensions. 

PACS Indices: 75.10.Jm, 75.40.Mg, 74.20.-z 
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The study of high-Tc superconductors^ continues attracting much attention. Re- 
cent progress has been made in the theoretical search for ground state superconduc- 
tivity in one band electronic models. The two dimensional t — J model has shown 
indications of d-wave superconductivity^ when analyzed near the phase separation 
regime^ at intermediate densities. It would be desirable to extend these observa- 
tions about superconductivity to more realistic three band versions of the Hubbard 
model, such as those introduced by Emery^ and Varma and collaborators.^ Not much 
is known for the superconducting properties of these models. Weak coupling and 
large N mean-field calculations applied to this extended model have indeed shown 
that a region of superconductivity exists near phase separation as in the one band 
case.^ However, these techniques are only approximate and should be supplemented 
by studies of the multiband Hubbard model with unbiased computational methods, 
like exact diagonalization and quantum Monte Carlo^'^. 

Varma and collaborators^ have suggested that the inclusion of a short-range 
Coulomb repulsion induces a charge transfer process, leading to the formation of 
tightly bound hole pairs on the oxygen ions. If these pairs are mobile they naturally 
lead to a superfluid condensate. Unfortunately, the charge transfer mechanism seems 
systematically correlated with a phase separation process, as was recently discussed 
in the atomic limit. ^ Phase separation occurs between a phase with a density of one 
hole per unit cell on the copper ions and a second phase which has a density of two 
holes per unit cell, with the charge located on double occupied oxygens. These pairs 
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are not mobile, and thus the system cannot become superfluid. In this sense, phase 
separation is an unwelcomed effect in this model. In spite of this problem, Sudb0 
et al.^° recently found indications that the one dimensional model exhibits supercon- 
ducting correlations immediately before phase separation. This conclusion is based 
on a study of the conformal field theory parameter Kp] Kp > 1 indicates that super- 
conducting power-law correlations are dominant in the ground state. (Unfortunately 
the actual pairing correlations in the ground state decay rapidly with distance^^.) 
While this result is very encouraging, it is difficult to predict in advance whether a 
region of superconductivity would exist near phase separation in the more realistic 
two dimensional problem. Thus, it would be desirable to have a model wich exhibites 
robust superconducting correlations in the ground state, and whose main features 
can be understood intuitively, allowing its extension to two dimensions where explicit 
numerical simulations are difficult. 

In this paper, a simple modification of the standard three band model is presented 
that addresses these issues. Even in the atomic limit, a clear separation between the 
parameter values associated with charge transfer and phase separation is evident. 
Including a hopping term, robust numerical evidence of superconducting correlations 
is observed. This three band Hubbard model for the Cu-0 cuprates is defined by,^'^ 

H = -tpd J2 i^laCja + c]^Cia) + UdY^ rii^riii + C/p X] ^i^'^jl + H + (1) 

<ij>a i j i 3 

where the operators are hole number and hole creation and annihilation operators, 
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Ud and Up are the Coulomb on-site repulsion energies for the copper and oxygen sites 
respectively, and and tp are the ion energies. The label i corresponds to Cu sites, 
and j denotes O sites and they are connected by a hopping amplitude tpd- The charge 
transfer energy is defined by A = ep — e^. It is generally accepted that Ud > A, Up 
in the cuprates, so we will work in this regime. The doping fraction is x = Uh/N, 
where Uh = (Nh — N) is the number of doped holes, Nh is the total number of holes 
and N is the number of Cu-0 cells. At half-filling, Nh = N, as in the insulating Cu-0 
compounds. Most of the results presented here correspond to low and intermediate 
doping {x < 0.66). 

As explained before, Varma et al.^ argued that the Hamiltonian Eq.(l) needs to 
be supplemented by a nearest-neighbor Coulomb repulsion, 

Hnn J2 '^i'^v (2) 

<ij> 

(where V is a positive number, and hole number operators) to induce the 

charge transfer mechanism, but this term also induces phase separation, at least in 
the atomic limit. To avoid this problem we have extended the range of the density- 
density correlations of Eq.(2). The relevance of the long range part of the Coulomb 
interaction between the holes, which prohibites phase separation was emphasized in 
a recent work by Emery and Kivelson.^^ In this work, experimental evidence was 
provided in support of the "frustrated phase separation" scenario which results in 
the presence of long range Coulomb interaction.^^ In particular, if the interactions 
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are of infinite range, then phase separation is certainly ehminated.^^ This leads us to 
replace Eq.(2) by a more general interaction 

= ^ 2. -1— e ^ , (3) 

where dim = — fml, and A is the range of the screened Coulomb repulsion. If 
A = oo, the interaction decays as the inverse of the distance, and it is of infinite 
range. Most of the results presented are for this limit. We also consider finite values 
of A to show that our conclusions are qualitatively valid even for finite range forces, 
provided that they operate over distances larger than a few lattice spacings. 

To motivate the introduction of a long-range term in the Hamiltonian it is conve- 
nient to first consider the atomic limit tpd — 0. Let us consider the one dimensional 
case first. In this case, it is not difficult to guess which are the possible states with 
lowest energy for different values of V. We have then checked that some of these states 
are effectively the ground state in a certain region of V using a simulated annealing 
algorithm described below. In Fig. la we show the ground state in the regime of small 
V. This state (denoted by I), has a hole at every copper ion, and the remaining of the 
holes are on singly occupied oxygens, distributed in a regular pattern to minimize the 
newly introduced long-range potential energy. In this region, for the two dimensional 
(2D) model, and for not too high doping, some antiferromagnetic (AF) correlations 
are expected to develope. On increasing V a phase separated state (denoted by II 
and shown in Fig. lb) becomes energetically competitive with state I.^ It contains a 
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large region of doubly occupied oxygens and another separated region with one hole 
per copper ion. These two densities coexist and changing the number of holes just 
increases or reduces their relative size. However, there is a third possibility which 
becomes energetically favorable due to the long-range Coulomb interaction. In Fig.lc 
we show a state (III) where all the holes are on doubly occupied oxygens. The reason 
for considering this state is that charge tends to spread uniformly over the lattice in 
the presence of long-range forces, rather than separating into different densities, as 
in the phase separated state. State III is a configuration with "preformed" pairs on 
the oxygens, and if this state becomes stable in the atomic limit, a finite hopping 
amplitude may render it superfiuid. In other words, once the pairs are preformed at 
zero hopping, it is reasonable to expect that they will behave as hard-core bosons 
when they are allowed to move, at least for a small hopping amplitude. Similar be- 
havior is observed in the attractive Hubbard model at large |U|/t. Note that state 
HI resembles a Wigner crystal of charge 2e pairs at the oxygens which is precisely 
expected to be stable due to long-range interactions. We have compared the energies 
of states I, II and HI numerically in the atomic limit for 12, 24 and 48 CuO cells. 
We have considered Ud = 7, Up = 1, A = 1.5, in arbitrary units, and x = 0.33. 
For the 12 cell lattice there is a small interval of V, [3.1, 3.7], where state II has the 
lowest energy. For the 24 cell lattice this interval is smaller and finally, for the 48 
cell lattice there is a direct crossing from state I to state III at ~ 3.38. Similar 
arguments are obviously valid in the two dimensional case, so our results are not 
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restricted to chains only. For illustration we have carried out an explicit calculation 
for a 4 X 4 Cu02 periodic cluster with four doped holes {x = 0.25). In this case 
the ground state is not obvious, and it is difficult to determine it by calculating the 
energy of each possible configurations because their number grow exponentially with 
the number of atoms. For this reason we developed a simulated annealing^^ program 
to find the ground state. This algorithm can be described as a sequence of Markov 
chains, each one generated at a fixed value of the temperature. The temperature is 
decreased between subsequent Markov chains. At each step of the Markov chain, the 
algorithm attempts to transform the current state into an state obtained by moving 
one hole to a nearest-neighbor site. The variation of the energy is computed, and the 
local change is accepted according to the Metropolis criterion. Using this technique 
we have determined the ground state of this cluster for the same set of parameters 
Ud, Up and A as before. For small V, the ground state has one hole per Cu atom 
and the doped holes are in O sites locates as far as possible from each other. As V 
is increased, a charge transfer from Cu to O sites starts to develope. For V ^ 3, 
a state with doubly occupied O sites separated from a region of single-occupied Cu 
sites becomes stable. For larger values of V, the ground state has the form illustrated 
in Fig. Id. These results in two dimensions give further support to our intuitive con- 
clusion that in the presence of strong long-range interactions, there is a phase with 
preformed pairs without phase separation. 

In order to find the A dependence of our results, numerical calculations where 
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carried out using the simulated annealing technique on chains of 12 and 24 CuO 
cells {x = 0.33). We have considered the same parameters as above. As an order 
parameter we used the number of doubly occupied oxygen sites. The results are 
shown in Fig.2a. We observed that our intuitive ideas about states I, II, and III are 
qualitatively correct, and they are dominant in a large region of parameter space. 
Only in narrow regions near the phase boundaries, especially for large values of A, 
did we observe other types of states become competitive. However, these states are 
simple variations of states I, II, and III (for example, including small fluctuations), 
so the essential features shown in Fig. 2 a are correct. Note that, for the 12 cell chain, 
there is a particular value of A ~ 8 where the phase separated regime becomes very 
narrow, and it may dissappear when the hopping parameter is switched on. (State 
III will reduce its energy by providing mobility to the pairs, but state II is too rigid 
to take advantage of this possibihty.) This value of lambda reduces to ~ 4 for the 24 
cell chain. We conjecture therefore that the phase separated regime does not exist 
for A > 4 or perhaps for even smaller values of A, so it is not necessary to have a 
fully unscreened 1/r interaction to obtain the effects described here. This value of 
A corresponds to two Cu-Cu lattice spacings (about 7.6^4 in the cuprates), and thus 
it is not physically unrealistic. The main features in Fig. 2 a are also present for the 
case of 8 doped holes {x — 0.66). Similar calculations for the 4x4 and 6x6 Cu02 
clusters, for x ~ 0.25, give essentially the same phase diagram, as it can be seen in 
Fig.2b. 
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The results summarized in Figs.l and 2a,b were obtained in the atomic hmit. 
Although they are very suggestive, the ultimate test of these ideas requires a finite 
kinetic energy term to give mobility to the carriers. The actual boundaries in pa- 
rameter space between different ground states will be determined by the competition 
between the Coulomb potential energy and the kinetic energy. To study this inter- 
play, and to establish the assumed existence of a superconducting regime, Lanczos 
diagonalization techniques were used to obtain the ground state of the Hamiltonian 
Eq.(l,3) explicitly on a finite chain. Due to the large number of states per unit cell 
the lattices that can be studied numerically are limited to a small number of cells. In 
this paper we studied 6 cells (12 atoms) with various hole numbers. For the results 
shown below we chose — 7, Up — 1 and A = 1.5 in units of tpd- These are rea- 
sonably realistic values for the cuprates.^^ We have also observed that our results are 
qualitatively similar over a broad range of parameters. 

First, let us numerically consider the issue of charge-transfer versus phase sep- 
aration numerically. A quantitative measure of charge transfer from Cu to O sites 
is given by the expected occupation of oxygen sites. However, a sharper indica- 
tion of charge transfer is given by the susceptibility associated with that quantity, 
XCT —< nQ > — < no >^, {no is the number of holes on the oxygens).^ The results 
for this susceptibility are shown in Fig. 2c for rih — 2 and = 4 on the 12 atom 
chain and a long range interaction. In both cases a peak is observed at a particular 
value of y 4. This peak presumably signals the onset of a charge transfer process 
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in this model. (Note also that xct in the region V > Ud shows additional structure, 
although this is an unphysical region.) xct however does not distinguish between the 
states II and III in Fig.l. To confirm that phase separation does not occur in the pres- 
ence of long-range interactions, we have measured the short wavelength component 
of the susceptibihty associated with the correlations of pairs of holes on O sites. 
This quantity (denoted by X) is normahzed to 1 for a fully phase separated state. In 
Fig. 2d, X is shown for the case n/j = 2 on the 12 atom chain. With a short range 
interaction, X —> 1 as V is increased. On the other hand, in the presence of long- 
range interactions, X goes through a maximum, and then converges to zero at large 
V, showing the absence of phase separation in this limit. The peak at intermediate 
V is due to the proximity of the phase separated state in the spectrum (as discussed 
previously in the atomic limit). At this point it is important to remark that the 
critical value of V (Vc) at which the charge transfer process occurs (namely ~ 3tpd) is 
physically acceptable. Some calculations suggest that tpd ~ 1.3eV in the cuprates,^^ 
and thus Vc obtained here is approximately 4eV. This is not far from estimations^^ 
of V based on measurements of the dielectric constant that suggest V ~ 2 — 3eV. 

Now let us investigate whether state III of Fig.lc becomes superfluid when the 
holes acquire mobility in the presence of LR interactions. For this purpose, we con- 
sider the pairing correlation, 

^M = ^E<^Vm^>' (4) 
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where the pairing operator is defined as = Cj^Cji, ^-nd j, j + m indicates O sites. In 
Fig.3a-b, the pairing correlations are plotted as a function of distance for two different 
densities, for tp^ — l,Ud — 7, Up — 1, and A = 1.5, and several values of V. For Uh = 2 
the correlations are robust at distance one which is the distance between two oxygen 
sites (~ 3.8 A in the cuprates). At larger distances the correlations decay rather 
rapidly, so our hypothesis is not clearly supported by these results (note however 
that these correlations are already stronger than those found in one band models 
near half-fiUing). On the other hand, with 16 holes on the chain (uh — 4) the pairing 
correlation is very strong even at the largest distances available on our finite clusters. 
The signal monotonically increases with V, for V < Ua-^^ 

The results presented thus far support a scenario in which s-wave (local) su- 
perconductivity appears with increasing V if LR interactions are included. To give 
further support to these claims, we have studied the dependence of the ground state 
energy on an external magnetic fiux To analyze this response, a phase factor 
gj*/iv -g introduced in the hopping term in Eq.(l).^'^° This is equivalent to allowing 
a magnetic fiux through the Cu-0 ring (chain with periodic boundary conditions). 
In Fig. 4, the ground state energy measured with respect to the energy at zero fiux, 
AE — £■($) — = 0), is shown as a function of $ for = 4 and a pure LR in- 
teraction (A = oo). For small V the energy has a single minimum at $ = (moo?27r). 
For V larger than a critical value the energy develops a second minimum at $ = tt, 
indicating the presence of mobile carriers with charge 2e in the ground state.^° This 
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anomalous flux quantization is in agreement with results obtained from the pair cor- 
relations. 

In summary, we have considered an extension of the standard Cu-0 electronic 
Hubbard-like model for the cuprates which incorporates long-range interactions. We 
found that in the atomic limit a charge-transfer process exists without phase separa- 
tion. The new dominant state in this regime has preformed pairs on the oxygen sites. 
With a hopping term that provides mobility to the holes, this state appears to become 
superfluid. The effects observed here are intuitively evident, and are valid both in 
one and two dimensions. These effects occur not only with a pure 1/r interaction but 
also for a finite-range Coulomb repulsion (A ~ 4) which might correspond to the real 
cuprate compounds. 
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Figure Captions 



FIG. 1 a) Ground state of the Hamiltonian Eq.(l) with long-range interactions for 
small V, in the atomic limit. The copper ions are represented by large circles 
and the oxygens by small circles. In this and the following figures in the atomic 
limit the alignment of spins is arbitrary. This state is ground state in the 
large 17^ regime if Up is nonzero, otherwise it becomes energetically favorable 
to move the holes of singly occupied oxygens into doubly occupied oxygens, b) 
Ground state of the model in the region of small A (short-range interactions) and 
large V. This state has phase separation and corresponds to the state previously 
discussed in Ref.[9]; c) New state discussed in this paper. It becomes the ground 
state of the problem in the case of large A (i.e. with long-range interactions) and 
sufficiently large V. This state has preformed pairs, and is not phase separated, 
d) Generahzation of the state shown in Fig.lc to a two dimensional lattice. 

FIG. 2 a) Phase diagram of the Hamiltonian Eqs.(l,3) in the atomic limit, obtained 
from a chain with 12 CuO cells (solid lines) and from a chain with 24 CuO cells 
(dotted lines), Ud = 7, Up = 1, and A = 1.5 (in arbitrary units), as a function 
of A and V. The doping fraction is x = 0.33. The boundaries were obtained 
numerically using a simulated annealing algorithm. "PM" denotes a region 
dominated by the states hke the one shown in Fig. la; "Ph. Sep." denotes the 
region where state Fig. lb dominates; and "Pairs" is the region with preformed 
pairs in the ground state, as exemplified by Fig.lc and Id. b) Same as a) for 
the 4x4 (solid lines) and 6x6 (dotted lines) and x = 0.25. c) Susceptibility 
XcT (see text) as a function of V for a 6 cell chain (12 atoms), tpd = 1 and other 
parameters as in Fig. 2a. The full squares denote Uh — 2 and the circles Uh — A. 
d) Order parameter X (see text) used to study phase separation versus V. The 
full squares denote A = oo, and the open circles are for the nearest- neighbor 
interaction. The number of holes is 8 {uh = 2). 

FIG. 3 a) Superconducting correlations C(m) (as defined in the text) for the long- 
range model at several values of V, on a chain with 12 atoms (6 cells) and 8 
holes (uh = 2). Note that m measures 0-0 distances. The full triangles denote 
results for V = 2, the open squares V = 4, and the full squares V = 6; b) As a) 
but with 10 holes (uh = 4). 

FIG. 4 The ground state energy AE{^) (measured with respect to the energy at zero 
flux) versus flux $ for a 6 cell-ring for various V and a long-range interaction 
(A = oo). 
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